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PI,  Changyong  Qin,  Benedict  College 
Co-PI,  Kevin  Huang,  University  of  South  Carolina 


The  following  topics  were  explored  by  experimental  and  theoretical  methods  in  the  funded 
period  of  08/2011-07/2014.  A  full  progress  report  is  presented  here. 


I  Interaction  of  oxygen  with  molten  carbonate 

1.1  DFT  study  of  C042-  and  C052'  in  gas  phase 

Using  DFT/B3LYP  methods,  two  intermediate  structures,  C042'  and  C052',  were  identified  and  analyzed. 
The  relevance  to  oxygen  reduction  at  the  presence  of  molten  carbonate  in  the  cathode  of  fuel  cells  was 
discussed.  The  binding  energy  of  C02  to  022"  is  calculated  to  be  143  kcal/mol.  The  high  binding  energy 
indicates  that  peroxide  plays  an  important  role  in  the  adsorption  of  C02  and  this  is  a  key  step  in  the 
reduction  of  oxygen  in  MCFCs.  The  formation  of  C052'from  C032'  and  02  is  exothermic,  with  an  energy 
change  of  -24.3  kcal/mol,  which  is  very  comparable  to  the  adsorption  energy  of  oxygen  on  the  surface 
defect  sites  of  the  LSM  cathode  of  SOFCs.  For  the  reaction  of  2C032'+02— >2C042',  the  equilibrium 
constant  at  873K  is  estimated  to  be  1.1  and  C042'  is  thermodynamically  favored  intermediate  in  the 
oxygen  reduction  process.  Through  orbital  analysis,  it  is  found  that  the  formation  of  C042'  will 
significantly  decrease  the  reduction  potential  of  oxygen  and  should  promote  the  ORR  process  in  SOFCs. 

More  details  of  this  part  can  be  found  in  the  published  paper  (Journal  of  Theoretical  and  Computational 
Chemistry).  Because  of  the  important  role  of  C042'  in  ORR,  next  step  will  be  concentrate  on  its  diffusion 
in  molten  carbonate  and  adsorption  to  the  LSM  cathode  surface.  On  the  cathode  surface,  C042'  will  then 
be  reduced  back  to  C032'  and  O2".  Finally,  O2'  will  be  transferred  to  electrolyte  and  the  ORR  process  is 
completed  in  the  cathode  of  SOFC. 

1.2  C042'  in  MC  as  reaction  intermediate  in  oxygen  reduction 

In  the  structure  of  (X2C03)4  X=Li,  Na,  K,  one  carbonate  was  substituted  by  C042'  and  each  structure  was 
re-optimized  using  B3LYP/6-31G(d).  The  substitution  of  C042"  does  not  significantly  change  the  structure 
of  cluster.  The  X-O  distance  is  1.968  A,  2.305  A  and  2.667  A  for  Li,  Na  and  K,  respectively.  The 
Mulliken  charge  decreases  in  the  order  of  Li,  Na,  K,  which  reflects  the  ionization  energy  of  the  alkali 
metals. 


In  this  part,  we  are  more  interested  in  the  structure  of  C042'  in  MC.  Compared  with  single  C042' 
calculated  in  the  gas  phase,  its  structure  remains  almost  same.  The  0-C032'  distance  is  about  1.465  A, 
similar  to  the  0-0  distance  in  peroxide. 

1.3  Oxygen  reduction  through  C042'  in  MC 

The  binding  of  oxygen  to  carbonate  forms  C052',  the  optimized  structure  of  C052'  is  very  close  to  that  in 
gas  phase  reported  and  only  negligible  change  was  detected.  This  C052'  will  then  react  with  another 
carbonate  to  form  C042'  as  shown  in  the  equations  below. 

2-  2- 

co  +  o  CO 

3  2  5 

2-  2-  2- 

CO  +  CO  2C0 

5  3  4 

For  step  two,  a  transition  state  was  located  on  the  potential  energy  surface  (PES)  for  each  type  of  MC. 
IRC  calculations  (Figure  2)  show  the  energy  barriers  and  confirm  the  reaction  path.  More  details  will  be 
included  in  the  manuscript  under  preparation.  The  shift  of  one  oxygen  to  the  other  carbonate  causes  the 
dissociation  of  oxygen  attached  to  one  carbonate,  which  makes  two  independent  C042',  which  then  serves 
as  oxygen  carrier. 


Figure  2:  PES  of  CC>5  +  CC>3  <->  2CC>4  (Top:  Li2C03,  Middle:  Na2C03,  Bottom:  K2C03) 

The  energy  barrier  of  oxygen  dissociation  to  form  C042"  is  estimated  to  be  49,  28,  41  kcal/mol  in  Li,  Na, 
K  MC,  respectively.  This  finding  needs  confirmation  form  experiments  and  also  further  modeling  with 
different  methods  will  be  performed  to  confirm  if  sodium  carbonate  will  benefit  the  ORR  process  more 
than  lithium  and  potassium  salts. 

1.4  New  Oxygen  Reduction  Mechanism  by  DFT  Modeling 

Combining  the  experimental  and  computational  results,  the  evidences  point  to  that  oxygen  reduction  is 
through  the  formation  of  C042'  in  MC.  This  active  oxygen  form  will  accept  electrons  and  release  an  oxide 


and  carbonate.  Further  studies  show  the  oxygen  can  easily  shift  between  two  neighboring  carbonates, 
which  means  the  oxygen  diffusion  in  MC  is  very  effective.  The  full  reaction  path  is  displayed  in  Figure  3. 
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Figure  3:  Formation  of  C042'  by  Oxygen  in  MC 


II  Fabrication  and  Test  of  SOFCs  with  MC  Modified  LSM  Cathode 

2.1  Fabrication 

Commercial  LSM  powder  (Nextech)  was  used  as  the  cathode  material.  The  powder  was  made  into  a  paste 
by  grinding  with  organic  binder  (V-006,  Heraeus).  The  paste  was  screen  printed  on  both  sides  of 
Lao.8Sro.2Gao.g3Mgo.17  (LSGM)  pellets  and  fired  at  1100  °C  for  2  h.  Li-K  binary  eutectic  salts  were 
prepared  for  infiltration.  Li2C03  (99%,  Alfa  Aesar)  and  K2C03  (99%,  Alfa  Aesar)  were  mixed  with  a 
molar  ratio  of  62:38  and  heated  at  650  °C  for  2  h  to  form  a  eutectic  melt.  In  addition,  eutectic  melts 
containing  0.5  mol  %  La203  (99.9%,  Alfa  Aesar)  or  0.5  mol  %  Gd203  (99.9%,  Alfa  Aesar)  were  also 
prepared.  For  infiltration,  the  eutectic  melt  was  ultrasonic  dispersed  in  ethanol.  A  few  drops  of  the  salt 
suspension  were  dripped  on  the  surface  of  LSM  cathodes  and  dried  at  100  °C.  This  step  was  repeated  till  5 
mg  salt  was  added.  The  samples  were  then  heated  to  650  °C  in  air  for  2  h. 

2.2  Electrochemical  test 

Silver  mesh  was  attached  to  the  LSM  cathode  using  silver  paste  (C8829,  Heraeus)  for  current  collecting. 
Electrochemical  Impedance  Spectroscopy  (EIS)  of  the  symmetrical  cells  was  obtained  with  an 
electrochemical  workstation  (IM6,  Zahner)  under  open  circuit  conditions  from  500  to  650  °C  with  50  °C 
interval.  The  frequency  range  was  swept  from  105  to  0.05  Hz  with  an  AC  amplitude  of  20  mV.  EIS 
spectra  were  simulated  and  fit  with  Thale  equivalent  circuit  software  (Zahner).  The  electrode  polarization 
was  taken  from  the  two  intersections  of  the  spectrum  at  the  highest  and  lowest  frequencies.  Thus  obtained 
polarization  resistance  was  then  divided  by  two  and  corrected  for  the  cathode  area  to  arrive  at  the  actual 
area-specific  resistance  of  the  cathode  polarization. 

2.3  Results  and  Discussion 

Figure  4  (a)-(d)  show  the  polarization  resistances  of  LSM  cathode  without  (baseline)  and  with  eutectic 
carbonate  melts  from  500  to  650  °C.  The  baseline  LSM  sample  exhibits  the  largest  polarization  resistance. 
After  infiltrated  with  5  mg  eutectic  melts,  polarization  resistance  has  been  decreased  significantly. 
Addition  of  rare  earth  oxides  La203  and  Gd203  further  improved  the  kinetics  of  oxygen  reduction  reaction. 
Among  all  the  tested  samples,  Gd203-added  Li-K  carbonate  melt  shows  the  lowest  area  specific  resistance 
(ASR),  a  57%  reduction  at  650  °C  and  80%  at  500  °C  compared  to  the  baseline. 


Figure  4  shows  the  Arrhenius  plot  of  cathode  polarization  for  base  and  eutectic  melt-infiltrated  LSM 
samples.  Introduction  of  eutectic  melts  into  the  porous  LSM  cathode  effectively  lower  the  activation 
energy,  suggesting  that  eutectic  carbonate  melts  can  promote  the  kinetics  of  oxygen  reduction  reaction 


Fig.5:  Temperature  dependence  of  area-specific 
polarization  resistance  for  the  baseline  and  eutectic 
melt-infiltrated  LSM  cathode  from  500  to  650  °C. 


From  the  reported  results,  we  can  conclude:  (1)  Molten  carbonate  can  improve  the  efficiency  of  LSM 
cathode  by  50%  in  SOFCs;  (2)  Addition  of  MC  to  LSM  can  reduce  the  activation  energy  of  oxygen 
reduction  in  the  cathode  of  SOFCs;  (3)  Rare  earth  oxides  La203  and  Gd203  can  also  improve  the  cathode 
performance.  However,  the  experiments  could  not  provide  any  details  on  the  fundamental  steps  relative  to 
the  enhancement  mechanism.  On  the  other  hand,  the  theoretical  team  already  proposed  an  oxygen 
reduction  mechanism  through  C042'  intermediate  and  this  is  consistent  with  the  results  from  experiments. 
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III  Fabrication  and  Test  of  SOFCs  with  MC  modified  LSCF/GDC  Cathode 
3.1  Fabrication  of  MC  modified  LSCF/GDC  cathode  cells 


A  commercial  LSCF+GDC  cathode  pink  (Sku:  232202,  Fuel  Cell  Materials)  was  used  to  screen-print  the 
two  identical  layers  of  porous  cathode  on  the  LSGM  electrolyte.  The  tri-layer  symmetrical  cells  were 
fired  at  1 100°C  for  1  h.  The  effective  surface  area  and  thickness  of  the  resultant  cathode  are  0.75  cm2  and 
30  pm,  respectively.  For  all  electrochemical  tests  below  650°C,  silver  paste  (C8829,  Heraeus)  and  silver 
mesh  were  used  as  the  current  collector. 

A  range  of  binary  Li-K  carbonate  compositions  was  prepared  by  melting  lithium  carbonate  (99%,  Alfa 
Aesar)  and  potassium  carbonate  (99%,  Alfa  Aesar)  at  650°C  for  2  hours.  The  resultant  melts  were  then 
broken-up  into  micron  size  powders  through  ball  milling,  which  were  subsequently  dispersed  in  ethanol 
under  ultrasonic  stirring  to  yield  a  suspension.  A  controlled  amount  of  the  suspension  was  then  uniformly 
distributed  over  the  cathode  surface  using  a  pipette,  followed  by  drying  at  100°C,  A  final  heat  treatment 
was  carried  out  at  650°C  for  2  hours  to  allow  carbonates  to  further  incorporate  into  porous  cathode 
microstructure.  The  compositions  of  the  binary  Li-K  carbonate  melts  was  selected  from  Li/K=  10/90  to 
90/10  (mol%),  including  the  eutectic  composition  at  62/38.  Influence  of  MC  loading  on  the  cathode 
performance  was  also  investigated  as  a  variable. 

The  Li-Na-K  carbonates  were  incorporated  into  in  the  cathode  scaffold  through  vacuum  infiltration  of 
constituent  carbonate  solutions  premixed  in  the  eutectic  composition  of  Li:Na:K=43. 5:3 1.5:25.0  (mol%). 
The  Li-Na-K  precursor  solution  had  a  concentration  of  lg/lOOml.  The  loading  of  carbonate  was 
controlled  by  the  numbers  of  infiltration.  After  each  infiltration,  the  samples  were  dried  at  100°C  in  air  for 
12h  to  evaporate  the  water  solvent.  Effect  of  MC  loading  on  the  cathode  performance  was  also  evaluated. 

3.2  Structural  characterization  and  cell  analysis 

The  spectra  of  Electrochemical  Impedance  Spectroscopy  (EIS)  was  collected  from  the  symmetrical  cells 
under  open-circuit  conditions,  and  further  analyzed  using  the  Thale  equivalent  circuit  software  (Zahner) 
to  extract  electrode  polarization  resistance  by  normalizing  it  to  the  actual  cathode  area  (0.75  cm2).  Typical 
impedance  scans  were  conducted  at  50°C  interval  from  650  to  550°C.  Twenty  minutes  were  allowed  to 
stabilize  at  each  temperature  before  the  scanning.  The  frequency  was  swept  from  105  to  0.01  Hz  with  20 
mV  of  stimulus  AC  amplitude.  A  long-term  stability  of  both  pristine  and  MC-added  cathodes  was 
simultaneously  evaluated  for  a  total  of  1,500  hours  at  600°C  with  a  side-by-side  cell  arrangement. 

The  microstructures  of  cathodes  in  all  the  symmetrical  impedance  cells,  either  pre-tested  or  post-tested, 
were  examined  with  a  field-emission  scanning  electron  microscopy  (FESEM,  Zeiss  Ultra)  equipped  with 
an  energy  dispersive  X-ray  spectroscopy  (EDS)  analyzer. 

The  electrochemical  performance  was  first  characterized  by  the  impedance  spectra  of  symmetrical  cells 
from  450°C  to  650°C  with  50°C  interval.  The  frequency  was  swept  from  105  to  0.01  Hz  with  20  mV  of 
stimulus  AC  amplitude.  Symmetrical  cells  with  varied  loading  were  compared  in  ambient  condition  to 
find  out  an  optimized  MC  amount  needed  in  cathode.  The  effect  of  oxygen  partial  pressure  on  the 
electrochemical  performance  was  evaluated  for  the  pristine  and  optimized  cathode  in  the  desired  ratios 
(latm,  0.21atm,  0.05atm  and  O.OOlatm).  A  flow  of  lOOsccm  was  passed  over  the  samples  at  all  the  time 
and  allowed  to  equilibrate  for  20min  in  each  atmosphere. 

The  same  cathode  was  applied  to  the  LSGM  electrolyte  supported  full  cell.  This  cell  has  an  anode  of 
Ni/GDC  cermet  with  thickness  of  50pm.  For  full  cell  test,  fuel  cell  was  sandwiched  between  alumina 
tubes  sealed  by  low  temperature  glass.  Voltage -current  and  impedance  measurements  were  conducted 
from  500  to  650°C  with  50°C  interval.  A  H2  flow  of  50sccm  and  an  air  flow  of  lOOsccm  were  supplied 
during  the  test. 

The  microstructure  of  cathode  from  both  pristine  and  MC  modified  cells  were  characterized  by  SEM 
before  and  after  electrochemical  measurement. 

3.3  Results  and  Discussion 
3.3.1  The  MC  loading  effect 


The  effect  of  MC  loading  into  the  microstructure  of  cathode  (LSM+GDC)  was  evaluated  using  EIS 
technique.  The  EIS  spectra  measured  in  the  temperature  range  of  500  to  650oC  are  shown  in  Fig.2.  The 
intercept  of  spectra  with  z’-axis  at  high  frequency  represents  ohmic  resistance  of  the  cell.  It  is  evident  that 
the  presence  of  MC  in  the  cathode  has  little  effect  on  it.  In  contrast,  MC  has  an  appreciable  effect  on  the 
polarization  resistance  of  the  cathode,  represented  by  the  length  on  the  z’-axis  between  the  high  and  low 
frequency  intercepts.  The  results  suggest  that  1.3  mg/cm2  is  the  optimal  loading  for  MC  to  yield  the 
lowest  polarization  resistance. 


Figure  6:  EIS  spectra  illustrating  the  effect  of  MC  loading 


3.2.2  The  effect  of  partial  pressure  of  oxygen  (p02) 

To  understand  the  charge  transfer  mechanism  of  MC  enhanced  ORR,  we  performed  EIS  measurements  as 
a  function  of  p02.  Fig. 7  shows  the  EIS  spectra  measured  as  a  function  of  temperature  and  p02  from  a 
pristine  sample.  As  is  evident  from  these  plots,  the  polarization  resistance  at  a  given  temperature  increases 
with  decreasing  p02;  the  ohmic  resistance  is  little  affected  by  the  p02  variations.  This  trending  suggests  a 
p- type  conducting  behavior  of  the  cathode.  In  addition,  it  seems  that  the  polarization  resistance  is  very 
sensitive  to  the  change  of  p02  for  this  type  of  pristine  cathode.  On  the  other  hand,  the  polarization 
resistance  as  well  as  ohmic  resistance  follows  typical  Arrhenius  relationship,  indicating  both  processes  are 
thermally  activated. 


Figure  7:  The  ESI  spectra  showing  the  effect  of  p02  on  the  pristine  cathode 

For  the  MC-modified  cathode,  the  polarization  resistance  shown  in  Fig.4  is  not  as  sensitive  to  p02  as  that 
of  the  pristine  sample.  The  shapes  of  EIS  spectra  are  also  different  from  those  of  the  pristine  sample. 
These  trends  suggest  that  the  charge -transfer  mechanism  in  the  MC-modified  cathode  is  different  from 
that  in  the  pristine  cathode. 

To  determine  this  difference,  we  need  to  construct  an  equivalent  circuit  to  simulate  the  EIS  spectra 
measured,  from  which  every  part  of  the  EIS  spectra  is  assigned  with  an  electrochemical  process.  The 
results  of  the  simulation  can  produce  a  series  of  parameters  for  the  electrochemical  processes.  From  these 
parameters,  the  charge -transfer  mechanism  can  be  interpreted.  This  part  of  work  is  currently  in  progress 
in  our  lab. 


Figure  8:  The  ESI  spectra  showing  the  effect  of  p02  on  the  MC-modified  cathode 


3.2.3  Microstructural  features 


The  microstructures  of  the  pristine  and  MC -modified  cathodes  were  also  investigated.  Fig.  9  shows  the 
microstructures  of  the  pristine  cathode  before  and  after  the  test.  There  is  no  difference  discernible  in  the 
SEM  images,  suggesting  that  the  microstructure  of  the  pristine  cathode  has  not  been  affected  by  the  high 
temperature  exposure. 

The  microstructures  of  the  MC-modified  cathodes  shown  in  Fig.  6  appear  to  be  affected  by  the  thermal 
exposure.  Before  testing,  the  molten  morphology  of  MC  is  clearly  seen  in  the  cathode  microstructure. 
After  the  test,  the  MC  phase  seems  to  disperse  into  much  finer  particles,  spreading  over  the  ORR  active 
sites.  The  presence  of  MC  at  these  active  sites  can  promote  charge  transfer  via  oxygen  absorber  (C042  ) 
and  transporter  (C052  ). 


Figure  9:  The  SEM  images  of  the  pristine  cathode  before  and  after  tests 
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Figure  10:  The  SEM  images  of  the  MC-modified  cathode  before  and  after  tests 
3.4  Fuel  cell  performance 

The  performances  measured  from  500  to  650°C  for  SOFC  with  and  without  MC-modified  cathodes  are 
shown  in  Fig.l  1.  It  is  evident  that  the  presence  of  MC  has  a  sizable  beneficial  effect  on  power  density.  At 
higher  temperature  such  as  650°C,  the  MC-modified  SOFC  shows  a  strange  behavior  at  high  current 
density.  We  will  investigate  the  cause  behind  this  phenomenon. 


Figure  11:  Comparison  of  V-J  curves  of  the  pristine  and  MC-modified  SOFCs 
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The  EIS  spectra  of  the  SOFC  shown  in  Fig.  8  further  confirm  the  V-I  results.  The  total  cell  resistance  of 
the  SOFC  with  MC -modified  cathode  clear  shows  a  much  lower  value  than  the  pristine  cell. 


Figure  12:  The  ESI  spectra  of  SOFCs  with  and  without  MC  modification 


IV  Raman  and  DFT  study  of  pyrocarbonate 

4.1  In-situ  Raman  spectroscopic  investigation 

The  Fi2C03  and  Na2C03  eutectic  mixture  (52:48  in  mol%)  was  first  synthesized  in  air  at  650°C.  After  a  2- 
h  hold,  the  melt  was  then  quenched  to  room  temperature,  followed  by  breaking  it  into  fine  particles  by 
ball  milling.  Thus  prepared  powders  were  then  packed  into  a  gold  crucible  that  was  subsequently  loaded 
into  a  high  temperature  stage  (Finkam  TS1500,  0°C~1500°C).  The  temperature  and  gas  were  controlled 
by  a  system  controller  (Finkam  PE95).  The  Raman  spectra  were  recorded  with  a  FabRam/HR  confocal 
Raman  system  (FabRam  Invers,  Horiba  Jobin-Yvon)  with  a  He-Ne  laser  operated  at  632.8  nm.  Since  the 
position  of  the  thermal  couples  in  the  high  temperature  stage  is  located  outside  the  crucible,  the  actual  and 
the  controlled  temperatures  of  the  MC  are  different.  The  melting  point  of  the  of  (Fi/Na)2C03  at  490°C  was 
used  to  calibrate  the  actual  temperature.  The  scattering  Raman  spectra  were  collected  in-situ  from  the 
carbonate  as  a  function  of  temperature  (in  the  range  of  RT-  600°C)  and  atmosphere  (in  N2,  air,  C02). 

4.2  DFT  modeling 

DFT  calculations  were  performed  at  the  B3FYP/6-31G(d)  level  using  Gaussian09  suite  of  quantum 
programs.  The  geometry  of  Fi2C205  and  Na2C205  were  optimized  first,  and  the  vibrational  frequencies 
were  then  obtained  from  analytic  second  derivatives  using  a  harmonic  oscillator  model.  In  addition, 
Raman  intensities  were  computed  by  numerical  differentiation  of  dipole  derivatives  with  respect  to  the 
electric  field. 

4.3  Temperature  dependence  of  Raman  spectrum 

Fig.l  (a)  shows  the  Raman  spectra  collected  from  a  eutectic  Fi2C03-Na2C03  (52  mol%  Fi2C03- 
Na2C03)  melt  over  a  band  range  of  650-1,850  cm'1  as  a  function  of  temperatures  in  a  pure  C02 
atmosphere.  For  the  carbonate  in  solid  state,  the  Raman  spectra  are  seen  to  contain  four  basic  vibrational 
modes  relevant  to  C032'  ions.  The  two  bands  at  1,078  and  1,094  cm'1  are  assigned  to  Vi  of  symmetric 
stretching  vibrations  in  Fi2C03  and  Na2C03,  whereas  the  bands  at  792  cm'1  and  870  cm'1  are  assigned  to 
v2  of  out-of-plane  bending  vibrations.  For  the  isolated  C032'  that  has  a  D3h  symmetry,  v2  is  Raman 
inactive.  However,  it  is  likely  that  v2  becomes  Raman  active  for  (Fi/Na)2C03  due  to  the  distortion  of  the 
C032'  structure  imposed  by  the  cations.  The  observed  bands  at  704  and  728  cm'1  are  assigned  to  v4  of  in¬ 
plane  bending  vibrations.  This  mode  is  a  double  degeneration  for  the  distorted  C032'  induced  by  Fi+  and 
Na+.  As  the  C032'  group  becomes  distorted  from  its  regular  planer  symmetry,  this  mode  splits  into  two 
components.  The  bands  at  1,375  cm'1,  1,404  cm'1,  1,531  cm'1  and  1,563  cm'1  are  attributed  to  a  split  v3  of 
the  asymmetric  stretching  vibrations  caused  by  the  existence  of  Fi+  and  Na+  around  the  C032"  ions. 
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Figure  13  (a)  The  Raman  spectra  of  (Li/Na)2C03  in  the  region  of  650-1850  cm'1  as  a  function  of 
temperature  in  C02  atmosphere;  (b)  Magnified  view  in  the  region  of  1,000-1,150  cm'1 

The  shift  in  Raman  band  with  temperature  is  better  viewed  in  Fig.l  (b),  a  magnified  spectrum  showing 
the  region  of  the  major  vrbands  at  1,078  cm'1  and  1,094  cm'1.  As  the  temperature  increases,  the  vrbands 
for  Li2C03  and  Na2C03  shift  toward  lower  wavenumber  and  eventually  merge  into  one  broad  peak  at  the 
melting  temperature  of  490°C.  This  shift  is  a  direct  result  of  lowered  force -constant,  elongated  C-O  bond 
length  and  weakened  Li(Na)-C-0  bond  strength  by  increasing  temperature.  When  the  temperature  reaches 
the  melting  point,  the  overtone  of  the  out-of-plane  bending  mode  (2*v2)  appears  at  1,762  cm'1.  It  should 
be  noted  that  a  broad  but  small  peak  near  970  cm'1  only  observed  in  solid-state  carbonates  appears  not 
directly  related  to  CC^2",  identification  of  which  is  not  possible  for  this  study.  We  speculate  that  the 
impurities  in  the  sample  could  be  a  source  for  this  unknown  peak.  A  detailed  assignment  of  the  measured 
Raman  bands  at  different  temperatures  and  atmospheres  is  summarized  in  Table  1. 

The  most  distinguishable  features  of  Raman  spectra  in  Fig.l  are  observed  when  the  carbonate  is  in  a 
molten  state.  The  four  bands  corresponding  to  the  v3  mode  disappear  from  the  spectrum  while  two  new 
broad  bands  at  1,317  cm'1  and  1,582  cm'1  emerge  within  the  same  band  width.  A  natural  question  to  ask 
is:  are  these  newly  emerged  Raman  bands  associated  with  the  C2052'  species? 

Table  1  The  Raman  frequencies  measured  and  assigned  for  the  eutectic  (Li/Na)2C03  at  selected 
temperatures  and  atmospheres 


Atmospheres 

co2 

Air 

n2 

Temperatures 

(°C) 

RT 

350 

455 

490  (melted) 

490 

(melted) 

490 

(melted) 

Frequencies 

704 

706 

706 

707 

707 

707 

(cm'1) 
n4  (weak) 

728 

720 

715 

- 

- 

- 

792 

790 

789 

790 

790 

790 

u2  (weak) 

879 

879 

881 

882 

885 

885 

1,078 

1,077 

1,077 

1,070 

1,072 

1,072 

Di  (strong) 

1,094 

1,090 

1,087 

_ 

_ 

_ 

1,375 

1,385 

1,386 

_ 

1,391 

1,392 

1,404 

1,415 

1,416 

_ 

1,421 

1,421 

n3  (weak) 

1,531 

1,520 

1,520 

_ 

_ 

_ 

1,563 

1,555 

1,545 

_ 

_ 

_ 

New  peaks 

_ 

_ 

_ 

1,317 

_ 

_ 

(broad) 

1,582 

Overtone 

- 

- 

- 

1,762 

1,762 

1,762 

(2*n2) 

Unknown 

- 

- 

972 

- 

- 

- 

4.4  Atmosphere-dependence  of  Raman  spectrum 

To  answer  this  question,  we  first  measured  Raman  spectra  in  different  atmospheres.  According  to  the 
enabling  electrochemical  reaction  shown  in  reaction  (3),  the  formation  of  C2052'  requires  a  source  of  C02. 
Fig. 2  compares  the  Raman  spectra  measured  in  N2,  Air  and  pure  C02  atmospheres  at  490°C  where  the 
carbonate  is  in  a  molten  state.  The  bands  at  1,072  cm'1,  790  cm"1  and  885  cm'1,  707  cm"1,  1,391  cm'1  and 
1,421  cm'1,  1,762  cm"1  shown  in  Fig.2  (a)  recorded  in  N2  atmosphere  correspond  to  the  symmetric 
stretching  (vi),  out-of-plane  bending  (v2),  in-plane  bending  (v4),  asymmetric  stretching  (v3),  and  the 
overtone  of  the  out-of-plane  bending  mode  (2*v2)  vibrations,  respectively.  The  Raman  spectrum  in  air  is 
almost  identical  to  that  collected  in  N2.  However,  the  peaks  at  1,317  cm'1  and  1,582  cm'1  shown  in  Fig.14 
(c)  are  only  observable  in  the  C02  atmosphere. 


Figure  14:  Raman  Spectra  of  C02  in  Different  Atmosphere 


The  strong  C02-dependence  of  the  bands  at  1,317  cm'1  and  1,582  cm'1  provides  a  crucial  hint  for  the 
formation  of  C2052'  via  the  C02  chemisorption  reaction  C02+C032'=C2052'. 

4.5  DFT  Modeling 

The  theoretical  support  to  the  formation  of  C2052'  is  provided  by  the  DFT  calculations.  The  calculations 
indicate  that  high  Raman  activities  of  C2052'  are  within  a  band  width  of  1,200-1,600  cm'1.  Specifically, 
the  active  Raman  bands  are  predicted  at  1,366  cm'1,  1,531  cm'1  and  1,566  cm'1  for  Li2C205  and  1,345  cm' 
l,  1,547  cm"1  and  1,579  cm"1  for  Na2C205,  respectively.  The  recorded  spectrum  in  the  band  width  of 
1,200-1,650  cm"1  from  the  MC  within  490-525°C  in  C02  atmosphere  could,  therefore,  be  an  overlap  of 
these  characteristic  Raman  peaks  of  Li2C205  and  Na2C205  in  this  band  region.  To  deconvolute  the  two 
unique  broad  peaks  around  1,317  cm'1  and  1,582  cm"1,  we  used  Guassian-Lorentizian  function  with  the 
six  theoretical  Raman  frequencies  as  the  standards.  The  results  are  shown  in  Fig.3,  where  the  black  and 
red  lines  represent  the  measured  and  modeled  spectra,  respectively.  Also  shown  are  the  individual 
spectrum  calculated  for  the  pure  Li2C205  and  Na2C205,  represented  by  pink  and  blue  lines,  respectively.  It 
appears  that  the  modeled  spectrum  is  dominated  by  the  Na2C03;  only  one  peak  at  1,566  cm"1  is  visible  for 
the  Li2C03  while  the  other  two  are  too  weak  to  be  seen.  This  is  primarily  due  to  the  differences  in  size 
and  polarizability  of  the  Na+  and  Li+  ions.  Overall,  the  measured  Raman  spectrum  is  in  excellent 
agreement  with  the  DFT  calculations  if  the  temperature  effect  is  factored  in. 


Figure  15:  Deconvolution  of  the  two  broad  bands  observed  at  1,317  cm'1  and  1,582  cm'1;  Inset:  the 
measured  and  DFT-modeled  Raman  spectra  in  a  band  width  of  1,200  -  1,650  cm'1 

4.6  Conclusion 

In  summary,  we  demonstrate  the  first  experimental  evidence  for  the  existence  of  pyrocarbonate  C2052' 
species  in  a  eutectic  Li2C03-Na2C03  melt  exposed  to  C02  atmosphere  through  a  combined  “DFT 
modeling”  and  “Raman  Spectroscopy”  approach.  The  existence  of  C2052'  species  is  a  key  support  to  a 
new  bi-ionic  transport  model  established  to  elucidate  the  fundamentals  of  the  high-flux  C02  transport 
phenomenon  observed  in  the  superior  mixed  oxide-ion  and  carbonate-ion  conducting  C02  separation 
membranes.  The  broad  Raman  bands  centered  at  1,317  cm'1  and  1,582  cm"1  are  a  characteristic  of  C2052' 
species  in  molten  carbonates  exposed  C02  atmosphere.  The  measured  characteristic  Raman  frequencies 
of  C2052"  are  in  an  excellent  agreement  with  the  DFT-model  consisting  of  six  overlapping  individual 
theoretical  bands  calculated  from  Li2C205  and  Na2C205. 


